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I. INTRODUCTION 


1 . 1 This document is the second quarterly technical report under Con- 
tract No. NAS J.2-2110 for "A Study of Information Processing and Trade-Off 
Analysis for Earth Applications Satellites." Under terms of the contract, it 
takes the place of the October monthly technical report. The work that has 
been outlined m the September monthly report is described herein m greater 
detail, as is all work accomplished since that report. Contract expenditures 
have been approximately 50 percent of the funds available . The work accom- 
plished has been successful, consistent with this amount of expenditure. 

1.2 Accomplishments m five areas will be discussed m detail in 
this report. 

a. Overall Simulation Concept 

b Time-m-Visual-View Model 

c. Earth Mode l/Sce nano Tape 

d. Extension of Radiometer Sensor Characterization 

e. Spacecraft Data Processing Model. 

1.3 In addition to the subprograms described in the first quarterly 
report^/ a computer subprogram has been generated to determine a schedule of 
satellite visual time m view, and flow charts have been prepared for each of 
the processing function subprograms within the spacecraft data processing 
model. In the next quarter's effort, the driver program for the spacecraft data 
processing model will be completed, thus the analysis effort for that portion 

Steen and J E. Thomas, Jr., Information Processing Simulation for 
Earth Applications Satellites, First Quarterly Technical Report. ORI 
TR 561, 8 August 1969. 
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of the simulation will have been completed, and a ground data processor will 
then be developed. Effort m the next quarter will also be directed toward the 
generation of an abbreviated scenario tape for use in simulation debugging and 
toward generation of computer subprograms within the spacecraft data process- 
ing model. 
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II. OVERALL SIMULATION CONCEPT 


2. 1 The objective of this contract is to develop a computerized tool for 
carrying out research on the systems aspects'of complex information manage- 
ment systems , typified by conceptually advanced earth applications satellite 
operations. A Digital Simulation Model is to be developed, de-bugged, and 
applied to a typical set of trade-off analyses involving such factors as onboard 
storage, onboard processing capability, communication limitations, ground 
facility depld£me~nt, etc. The iritent of this contract is not to plan future 
missions (although useful data for such planning may well result from this 
contract) but rather to develop a versatile tool and technique for studying 
information systems of general type which may be required in future NASA 
missions. 

2.2 In order that it remain consistent with these user requirements for 
the simulation, the operational concept has been modified (from the single large 
simulation originally envisioned) to an '‘on-line" -"off-line" concept, m which the 
"on-line" portion of the simulation is used to satisfy user requirements by opera- 
ting with a scenario generated earlier using the "off-line" portion of the simula- 
tion. A gross flow diagram of the simulation organization and models is pre- 
sented m Figure 1. 

2.3 The work accomplished to date has been guided by the philosophy 
that a broad simplified version of the simulation should be initially developed, 
to be followed by increased sophistication m critical areas of the simulation 
determined by the experience gained m the early stages of the study. In addi- 
tion, an attempt has been made to configure the simulation such that the input 
data requirements are those data that— 

a. Are subject to frequent change as various 
trade-offs and design alternatives are being 
evaluated 
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Slide Rule 


'Off-line" Models 



FIGURE 1 . INFORMATION PROCESSING SIMUIATIOI^OPERATIONAL CONCEPT 












b. Gan be obtained or estimated without extraordinary 
effort on the part of the user. 


MODELS 


2.4 The following paragraphs discuss briefly each of the models utilized 
by the simulation. More detailed treatment of the models for which analysis 
has been completed are contained m later sections of this report. 

Time-in-Visual-View Model 

2.5 The time-m-visual-view computation utilizes the ground station 
coordinates (latitude, longitude) and the satellite ephemens characteristics 
(altitude, period, inclination angle) to compute a schedule of time m visual 
view for a specified number of hours, days, or orbits. This time -m-visual- 
view schedule will be subsequently used m the computations of the character- 
istics of the communications link. Analysis and programming of this model 
have been completed. Section III of this report presents a more comprehensive 
discussion of this simulation task. 

Down-Link Data Transmission Model 

2.6 The down-link data transmittal model consists of two routines, an 
RF routine and an error rate routine. The RF routine utilizes the time-m-view 
schedule which has been previously developed to compute the cam er-to- noise 
ratios for the spacecraft-to-ground-station communications link during the times 
indicated by the time-m-view schedule. The communications link character- 
istics required as input are those required to perform a normal path loss/power 
budget computation. The computation has been analyzed and programmed. 

2.7 The error rate model considers both the increase m data rate due to 
coding for error reduction purposes and the anticipated error rates as a function 
of the computed camer-to-noise ratios. As presently envisioned, the compu- 
tation of error rate versus camer-to-noise ratio would be accomplished by means 
of a table look-up m the program, with the entries m the table varying as a 
function of the error coding used. However, another method of computing error 
rate has been considered, it is not yet sufficiently well defined for presentation 
m this report. The time m view, bandwidth, and total data available for trans- 
mission are then used to determine the percent of data dumped during a pass 
over a ground station. 

Earth Model 

2.8 The earth model will provide synoptic detail for such items as cloud 
cover, land/water interfaces, icebergs, and agricultural areas. Because a con- 
siderable amount of flexibility is desirable in modeling earth conditions, the 
data will be manually produced through the use of a slide rule (a clear plastic 
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slide containing the subsatellite points for the orbital period and inclination 
angle of interest, superimposed upon a Mercator projection of the earth). 

Section IV of this report presents a more comprehensive description of the 
earth model. 

Scenario 

2.9 The scenario is a multi-dimensional table containing the earth and 
satellite/earth events of interest and keyed to the longitude of the ascending 
node. The longitude of the ascending node is partitioned into bands within 
the array. These bands can be of any desired width, thus providing the capa- 
bility of including a greater or lesser amount of detail m the earth model, at 
the penalty of increased storage requirements in the simulation. Section IV of 
this report presents a more comprehensive description of the scenario. 

Sensor Models 

2.10 The sensor models utilize sensory characteristics (camera sensors — 
focal length, image resolution, etc.; radar sensors — pulse width, beamwidth, 
power, etc.); and sensory duty cycles (time profile), as well as the satellite 
ephemens characteristics, m order to compute information relative to sensory 
ground coverage (coverage area, ground rate, digital resolution, bandwidth). 

Thus, the user will specify a complement of sensors and their characteristics, 
as well as a time profile of their "on-off" -cycles, and the simulation will com- 
pute the data rate and the total data gathered over the time period of interest. 
Characterizations have been accomplished and programmed for passive imaging 
sensors, passive non- imaging sensors, and two types of active imaging sensors. 
Detailed descriptions of sensor characteristics were contained m Section II and 
Appendix B of the First Quarterly Report. Extensions to these characterizations 
are found m Section V of this report. 

Spacecraft Data Processing Model 

2.11 The output of the sensory models feeds into the spacecraft data pro- 
cessing model. Figure 2 illustrates additional details of this model. As pre- 
sently envisioned, the user will specify, on input to the on-line portion of the 
simulation, the operations required to be- performed upon the output of each 
sensor. Additionally, he will specify the state-of-the-art characteristics (for 
the time period of interest) of the equipments to accomplish, these operations. 

The operations and the characteristics specified must be in correspondence to 
those which the simulation can recognize (see Figure 2), but it is not necessary 
to specify all of the available operations . 

2.12 The computations performed sequentially for each sensor within the 
spacecraft data processing model reflect the changes m the data rate due to 
the specified processing operations, the storage required due to these data 
rates, and the time delays that are introduced due to these processing operations. 
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FIGURE 2 


SPACECRAFT DATA PROCESSING MODEL 











2.13 Additionally, the data rates from each sensor, after each operation, 
are compared to the input capacity specified for the next operation. Thus, for 
example, if multiplexer operations are specified as having an input capacity of 
200 kilobit s/sec , and if Sensor 1 has a continuous information rate of 100 
kilobits/sec, then that sensor will utilize one-half the capacity of a single 
multiplexer unit. The model keeps track of the time relationships of the in- 
dividual sensors vis-a-vis one another and allocates the required number of 

• units for each operation based upon the total load of all the sensors. At the 
completion of each simulation run, the spacecraft data processing model speci- 
fies the total number and the connectivity of each type of unit required m order 
to perform the required operations specified m the simulation input. Additionally, 
the composite bit rate, storage requirements, and- time delays are computed. 

2.14 Although this concept represents only a minor shift m emphasis from 
the original concept (Proposal Subtask 1.3), the present concept has advantages 
of the original concept as follows: 

a . The output data rates of the sensors are not 
immediately obvious from the sensor specifications . 

Thus, the processing configuration specified a prion 
may not be adequate to process the sensor output data, 
and several modifications and reruns may be required 
before convergence. The present concept will provide 
a workable configuration with a single simulation run. 

b. The original concept required specification of processing 
configurations m a fair amount of hardware detail. For 
the simulation to be sensitive to these configurations, 

it would be necessary to know information such as the 
number of floating point adds for a particular data com- 
pression algorithm. It is considered that generation of 
this type of information would put a severe strain on the 
user of the simulation to the point where the simulation 
would receive limited use. The present concept operates 
on a functional basis rather than a hardware basis (data 
compression unit versus control unit, I/O, CPU, Mam 
Memory, etc.) and is thus able to operate at a higher 
indenture level with a resultant side-stepping of these 
"gritty details." It should be noted that at some point 
m the design process it will be necessary to develop 
the details of the software package for data processing 
algorithms, but these are beyond the scope of this 
contract . 
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Ground Station Data Processing Model 

2.15 Conceptual development of the ground station processing model is 

underway; however, the accomplishments m this area are not sufficiently 
developed to permit detailed discussion at this time. 
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III. TIME-IN-VIEW/ COMMUNICATIONS LINK MODELS 


TIME IN VIEW 

3.1 In Section IV of the first quarterly report, ^the analysis and result- 
ing maximum time-m-view model are described. The model produces a maxi- 
mum time m visual view, given satellite ephemens characteristics and earth 
station location information, and is very useful m establishing an upper bound 
for data transmission times between satellite and ground station. 

3.2 Such information is useful m the selection of appropriate orbits and 
the establishment of broad design criteria such as gross estimates of satellite 
storage requirements and data link characteristics . However, for scenario 
generation purposes it is necessary to produce actual, rather than maximum, 
time-m-visual-view data and, further, to produce a schedule of time-m- 
visual-view events . As formerly designed, the maximum-time -m-visual-view 
model was not capable of producing this information, and it was modified to 
fulfill these more stringent requirements . 

3.3 Analysis conducted as part of the modification effort indicated that a 
major portion of the maximum -time -m-visual-view. model was usable without 
modification if it could be supplied with a pomt(on the satellite orbit) within 
the earth station cone of visibility caused by earth station look -angle 
limitations . 

3.4 Subsequent analysis directed toward determining the .most efficient 
method of finding the appropriate point on the satellite orbit indicated that an 

■^P. Steen and J. Thomas, Jr., Information Processing Simulation for Earth 
Application Satellites, First Quarterly Report , ORI TR 561, 8 August 1969. 
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iterative rather than an analytical solution should be used, the rejection of an 
analytical solution was primarily a result of the limitation of spherical 
trigonometry, i.e., the multiplicity of solutions to trigonometric equations and 
the spherical trigonometric angular limitations tr). The iterative solution is 
a simple two-step procedure: (a) determine the satellite orbital point at which 
earth station and satellite longitude are coincident and (b) determine the 
closest satellite approach to the earth station coordinates. Both procedures 
use equations previously developed (Section IV of the first quarterly report)-^/ and 
make use of common iterative procedures to find the desired solution; there- 
fore, they do not require further detailed treatment. Appendix A contains a 
complete flow diagram for the modified time-in-visual-view routine. 

3.5 Figure 3 presents a sample output from the time-in-visual-view 

model The printout presents a time-m-visual-view schedule m time of arrival 
order, with one list for each of the two first ascending node longitudes con- 
sidered. The data represents a typical time-m-view schedule for approximately 
a 1-day period, i.e., 14 satellite orbits. 

COMMUNICATIONS LINK 


3.6 This section describes the communications link computations which 

have been implemented. These computations supply the slant range and the 
look angles used m the carrier- lo-no is e (C/N) ratio computation. Figure 4 
illustrates the angles and distances used m the' calculations link. Consider 
a satellite at an altitude h s and a ground station at the earth's surface. The 
satellite, ground station, and center of the earth define a plane triange with 
the angle fi at the vertex corresponding to the earth's center. This angle can 
be shown to be 


£ = cos 1 Ceos LAg cos LA Qg cos | LOg - LOqq | + sm LA^g sm LA 3 ] 


LAg = 

satellite latitude 

L° s = 

satellite longitude 

la gs 

ground station latitude 

lo gs = 

ground station longitude 


The distance between the ground station and the satellite is given by 
D = C(R + hg ) 2 + (R ) 2 - 2 (R+hg) (R) cos/3]® 
where R is the radius of the earth. 




Ibid. 
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1 FAIRBANKS , ALASKA* 

2 ROSMAN, NORTH CAROLINA* 

3 BARSTOW, CALIFORNIA 

4 CANBERRA, AUSTRALIA 


* PRIME STATIONS FOR ERTS -A AND -B 


FIGURE 3 . TIME-IN-VISUAL-VIEW SAMPLE OUTPUT - SAMPLE A 




EARTH STATIONS 

IDENT NO LOCATION 

1 FAIRBANKS, ALASKA* 

2 ROSMAN , NORTH CAROLINA* 

3 BARSTOW, CALIFORNIA* 

4 CANBERRA, AUSTRALIA 

* PRIME STATIONS FOR ERTS -A AND -B 


FIGURE 3 (Cont) - SAMPLE B 
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Next the angle a is computed from 


~1 rR sin £1 
a = sin l— g-j. 

This is the look angle from satellite to ground station, which is used to compute 
the satellite antenna gam as a function of the angle off the boresight line. The 
ground station antenna gain is computed for the angle (a + /S), which is the 
deviation from the local vertical at the ground station. Alternately for directable 
antennas the maximum antenna gam can be used. 

3.7 The carrier- to- noise power ratio at the receiver can be expressed as 

(transmitted power) (antenna gams) 

_ P carrier _ l° sses 

P noise noise power 


Transmitted power is supplied as a direct input to the program. The antenna 
gams are entered as a tabular look-up of gam versus angle off the boresight 
and are computed as a function of the angles previously discussed. Equip- 
ment losses are supplied as input to the program. Path loss is computed by 
the express ion^/ using the slant range previously computed. 

Path loss = (6.03) (10 s ) f 2 D 2 
where f = frequency (MHz) 

D = slant range (nmi) 


3 . 8 The noise power is computed as 

noise power = F • K • T • B 

where F = receiver noise figure (units) 

K= Boltzman's constant (1.38 x 10 j/°K) 

T = antenna temperature (°K) 

B = receiver noise bandwidth (MHz) 


3 9 The computed- C/N ratio-is next compared to a threshold value 

(specified as input data) to determine if the communications link is suitable 
for data transmission. The computation is done for both the up-link and the 
down-link, since their C/N ratios will generally be different. Figure 5 pre- 
sents a flow chart of the computation. This computation has been run separate- 
ly on the computer but has not as yet been integrated with the "Time-m~ 
Visual-View" program. 


3/ 


Radiation Inc , First Quarterly Report for Proiect IRIS, June 1965. 
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FIGURE 5. SATELI[XE/ ground station 
COMMUNICATION. LINK COMPUTATION 
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OPERATIONAL USE 


3.10 Having specified the required input parameters, the computer pro- 
grams will first compute a schedule of times in visual view (i.e. , line of 
sight) for each ground station location for as many orbits as desired. This 
computation is accomplished sequentially for each ground station location. 

3.11 The schedule of time m view and time out of view is next used to 
compute the spacecraft to ground station C/N ratio as discussed above 

3.12 It is anticipated that schedule(s) will be generated for particular 
configuration s) of ground station locations and communications link para- 
meters and that the output of these programs will be used as a .portion of the 
input for the scenario tape used to drive the "on-line" simulation. 

Input Data Required 

3.13 The following input data will be required for the time-m-view/ 
communications link model: 

a Location of all ground stations (latitude and longitude) 

b. Satellite ephemens data (assumed circular orbit) 

Height (nmi) 

Period (mm) 

Longitude of ascending node (deg) 

Inclination angle (deg) 

c . Communications link data 

Frequency (MHz) 

Transmitted power (W) 

Receiver noise bandwidth (kHz) 

Receiver noise figure (dB) 

Antenna temperature (°K) 

Antenna gams (dB) 

System losses (dB). 
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IV. EARTH MODEL/S 0ENARIO 


4.1 It has become apparent that there is considerable interest m model- 
ing earth conditions m a good deal more detail than was originally anticipated. 
Specifically, simulation of such items as cloud cover and land/w.ater interfaces 
is desirable m synoptic detail, because sensor on/off must be related to these 
events. In addition, the earth model must provide considerable flexibility m 
modeling earth conditions- because of the diversified utilization projected for 
the information: processing simulation Furthermore, since the scenario will be 
required to store the outputs of the earth model, it, too, must provide similar 
flexibility. 

SCENARIO 

4.2 The scenario will consist of a multi-dimensional table of earth 
events of interest, with entries m the table keyed to the longitude of the ascend- 
ing node. Figure 6 illustrates preliminary notions regarding the type of entries 
which would be included in the table. Consulting Figure 6, one sees that the 
longitude of the ascending node is partitioned into bands of longitude. These 
bands can be of any desired width, thus providing the capability of including 

a greater or lesser amount of detail m the earth ‘model at the penalty of increased 
storage requirements in the simulation. 

4.3 The entries m the table represent minutes into the orbit drawn from 
a 600-nmi orbit with a period of approximately 107 min. Thus, m Figure 6, it 
is specified that, for orbits which have an ascending node longitude between 
60°,E and 70°E, the subsatellite points will be over open water during minutes 
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Long, of 
Ascending Node 
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Water 

Land/Water 

Interface 

Icebergs 
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Field 

Clouds 

a * L< b 
c £ L< d 

600 £ s L< 70° E 

0-6, 23-26 5, 

6-3, 20-23, 

23—26 . S , 

38-42 

* 

' 

43-79, 86-101 

26.5, 30-33, 
43, 79-80, 86, 
101-104 


30-33 




*Place as desired for particular scenario^ 



FIGURE 6. SCENARIO TABLE 








0-6, 23-26. 5 , 43-79, and 86-101, additionally, wheat field investigations can 
be conducted during minutes 38-42, etc The event categories which have been 
specified are: 

a . Satellite over open water 

b. Satellite over land/water interfaces 

c. Satellite over iceberg areas 

d. Satellite over wheat field areas 

e. Satellite over cloudy areas. 

Additional events can be specified as desired. 

4.4 Given the existence of such a table m memory, the sensor on/off 
cycles can be specified on input to the program as keyed to the occurrence of 
particular events of interest. Thus, for example, if a camera system with 
adaptive resolution and frame rate is hypothesized, it will be possible to specify 
a high-resolution, rapid-frame-rate mode during events relating to land/water 
interface investigations, and a low-resolution, slow- frame-rate mode during 
events relating to open water investigations. 

4.5 It remains to determine a method of easily generating the event times 
necessary for the scenario. Figure 7 illustrates a method for generating these 
event times in a rapid, efficient manner. With reference to Figure 7, the sub- 
satellifce points for an orbit of approximately a 107 -mm period are superimposed 
upon a Mercator projection of the earth. The numerical entries in Figure 6 have 
been drawn from Figure 7. Thus, one sees that the entries for iceberg investi- 
gations correspond to subsatellite points over artic waters, and the times corre- 
sponding to wheat field investigations correspond to subsatellite points over the 
Western United States. 

4.6 Pursuing this thought further, there is a certain appeal (m terms of 
reduced computing time and ease of use) m adding two more events to those 
shown m Figure 6, i.e. , day/night and satellite/ground station communications 
link. If this were done, it would not be necessary to rerun the time-m-view 
and communications link models for investigations m which the number and 
location of ground stations are not changed. These models could be run for 
several complements of ground stations, and the results could be added to the 
scenario/scenarios. Day/night could be handled by a shading of portions of the 
earth model slide rule — or by a separate, appropriately shaded slide. 

4.7 One can then envision a few standard scenario tapes supplemented 
by a data deck for scenario generation m which a particular scenario is changed 
by selective modification to portions of the data deck cards. Thus, for example, 
if it is suspected that the results of a particular investigation are sensitive to 
the complement of ground stations used for that set of runs, a set of cards repre- 
senting times m view for a different complement of ground stations could be sub- 
stituted in the data deck for scenario generation. The sensory models and the 
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SUB POINT TRACK 



FIGURE 7 . SLIDE RULE ILLUSTRATION 







processing models would then be rerun against the new scenario tape and the 
results compared. Similar reasoning applies to cloud cover or, indeed, any 
entry m Figure 6. 

EARTH MODEL 

4.8 It is considered that a slide rule, with a transparent slide of the sub- 
satellite points, that could be slid along a Mercator projection of the earth 
would provide a means for rapdily developing scenarios with as much granularity 
as desired (20°, 10°, 5°,etc.). 

4.9 A slide rule of the type described, which can provide the desired 
flexibility, has been developed for use in the manual entry of earth events into 
the simulation m the detail required by the particular situation at hand. Such 

a slide rule has been designed for a 600-nmi orbit and is contained m Appendix 
B. Due to the reproduction process used for this report, the scales of the Mer- 
cator projection and the slide rule are not precisely identical, however, a more 
precise slide rule will be contained m the final report. 

4.10 A different slide would be required for each altitude and inclination 
angle of interest These could then be used to generate the scenarios which 
would then be transferred to magnetic tape. For the present, only one scenario 
will be developed. The amount of work required to integrate such a scheme into 
the simulation remains to be evaluated. Further development of the details of 
this approach should reveal any pitfalls . 
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V. EXTENSIONS OF SENSOR CHARACTERIZATIONS 


5.1 In. evaluating the radiometer characterization m the first quarterly 
report,— m light of current simulation objectives, it was apparent that ex- 
tension of the radiometer characterization to specify explicitly the relation- 
ships between the parameters m that characterization m terms of more detailed 
hardware parameters was indicated. In the following paragraphs, discussion 
of the work accomplished m this area is applicable to both infrared and micro- 
wave radiometers . Much of the extension analysis, has been drawn from 
Horan's work m the July 1968 issue of Spectrum.— 

5.2 The data rate for a passive radiometer has been characterized as a 
function of the radiometer scan rate, scan width, and resolution, as well as 
the satellite velocity and height.—' The analog data rate, which is simply the 
number of resolution elements, N, per unit time, is given m the following 
expression: 

N = — S = ■ S. (5.1) 

V fVT 


— P. Steen and J. Thomas, Information Processing Simulation for Earth 
Applications Satellites, First Quarterly Technical Report, ORI TR 561 , 

8 August 1969 . 

2 / 

— J. Horan, "Spacecraft Infrared Imaging," IEEE Spectrum , July 1968. 

3/ 

— University of Michigan, Peaceful Uses of Earth Observation Spacecraft, 
Vol. Ill, Sensor Requirements and Experiments, NASA, CR-588, 

3 October 1966. 
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where 


N is given m resolution elements/second 

a = scan width (sweep angle) 
y = angular resolution 
S = scan rate 

— instantaneous solid angular field of view . 


For no overlap, the scan rate S can be shown to be 



where v = satellite velocity 
h = satellite altitude* 


(5.2) 


With overlap, the scan rate is simply modified by the overlap coefficient, 
P, i.e . , 


S = 


v 

yhp 

i 


(5.2a) 


The effective scan angle is given by the following expression: 


data time + flyback time , , T d + 
e data time ' 7 


T r 


a 


and it is seen that 


(5.3) 


Thus 


data time + flyback time ~ 



(5.4) 



(5.5) 


where Tf = flyback time. Therefore, when the flyback time has been 
determined, the effective scan angle (<y e ) will be used for determining 
the data rates . 


Substituting m the expression for analog data rate, N, yields 



or for overlap conditions 



(5.6) 
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(5.6a) 



x 

where the angular resolution, y, = {0)°, the square root of the instantaneous 
solid angular field of view . 

The bandwidth, hX, required for this analog data rate is equal to 

Af = N/2 (5.7) 

Now, the radiometer signal = P det * 41 (5.8) 


where* P det ” P° wer focused on the detector 

Si — detector sensitivity (responsivity) m volts out/ 
watts m 

P det = A ° n77Nr (5,9) 

where Aq = area of collecting optics 

O = instantaneous solid angular field of view 
rj - optical efficiency 
Nt= radiance of target. 

(N) Radiometer noise = Johnson noise of detector resistance 

= J 4KTRAf ( 5 * 1 0). 

where K = Boltzman's Constant 

. T = temperature 

R = detector resistance 
Af = noise bandwidth. 

Signal/noise ratio (SNR) = ^ (5.11) 

/ 4KTR Af 

Therefore, SNR~ - - Q— ~ an d analog data rate (N) ~ , thus we have - 

(Af)e ZU 

established relationships which allow trade-offs between the radiometer 
data rate, which is a prime parameter of the simulation, and a wide variety 
of hardware parameters, i.e. , S/N ratio, bandwidth, field of view, area of 
collecting optics, etc. Further, it can be shown that—/ 

£ T= _EL_ (5 . 1.2) 

■J Bt 


4/ 


Ibid. 


25 


where 


AT — receiver minimum detectable temperature difference 
(thermal resolution) 

K = system constant = f (modulation scheme) 

T = system noise temperature 
B = equivalent-square predetection bandwidth 
T — post-detection integration time 

and for a Rayleigh limited system and S/N = 1 


where 


AT ph 


CD 
1 . 2 \ 


(¥ )* 


( 5 . 13 ) 


C = system figure of merit = 


time to scan a line 


total time from line to line 
D = diameter of collecting aperture 
X = wavelength of operation 
a = scanning angle 

v/h = velocity/height ratio of the space vehicle. 


Td 


T d +T f 


Thus, the data rate has been shown to be a fQ therefore, trade-offs can 

be made between thermal resolution and data rate m the simulation, also band- 
width and integration time. 



VI. SPACECRAFT DATA PROCESSING MODEL 


6.1 This section reports the progress made to date m the development of 
the spacecraft data processing model Work accomnlished includes 

the establishment of the model modus operandi as well as the development 
of detailed specifications for each processing function within the data 
processing model and the production of detailed flow charts for each pro- 
cessing function. Preliminary flow charts are contained in Appendix C. 
Additional effort has been directed toward detailed specification of the mini- 
mization procedure to be utilized within the model driver; however, results 
m this area are not sufficiently advanced for reporting at this time . 

MODEL ORGANIZATION 

6.2 As presently envisioned, the spacecraft data orocessing model 
will provide the user with seven processing functions from which the desiied 
processing configuration for each sensor may be selected. These are: 

a. Signal conditioning 

b. A/D conversion 

c. Multiplexing or commutation 

d. Buffering 

e. Overlap redundancy removal 

f. Information preserving data compression, and 

g. Entropy reducing data compression. 

6.3 The objective for the spacecraft data processing model is to con- 
figure a data processor, using the fewest possible number of processing 
functions required to process the stream of data emanating from the sensors, 
while not exceeding pre -specified state-of-the-art limitations for each 
processing function. To fulfill this objective, the model will be organized 
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into a driver routine and seven processing function routines, one for each 
processing function. 

6.4 With this objective m mind and recognizing that the spacecraft 
data processing model is designed for use as a planner's tool, state-of- 
the-art compliance of each processing function is of primary importance 
within the model . The rationale for actions to be taken relative to state-of- 
the-art compliance can be categorized into three basic cases: 

a. If the requested operation (requested by the driver 
routine) requires less than the state-of-the-art 
capability for a particular function, the processing 
function will be considered to have excess capacity 
to be used for additional processing. 

b. If the requested operation requires capability m 
excess of the state-of-the-art capability for a 
particular function, corrective action will be initiated 
to reduce the requested capability to within the state- 
of-the-art. 

c. If the requested operation cannot be reduced to within 
the state-of-the-art capability, the excess capability 
requested will be determined and supplied with the 
model output for further evaluation. 

6.5 While the processing required to accomplish the objectives 
of each of ihese three cases will be contained within the driver routine, 
since that routine has access to the a prion information required, the 
processing function routines will accomplish the state-of-the-art com- 
pliance testing, and, therefore, must supply the driver routine with the 
state-of-the-art compliance data. Flow charts for the seven processing 
functions are presented m Appendix C. 

PROCESSING FUNCTION SPECIFICATIONS 

6.6 The following paragraphs describe the rationale for each processing 
function, the assumptions made, and the methodology to be used m the 
model for each processing function. Table 1, which summarizes these data, 
is presented at the end of this section* Each processing function, however, 
requires (within the state-of-the-art limitations) the specification of a time 
delay which has the effect of shifting the input time profile by the amount of 
the delay. In addition, an output data transmission time is computed for 
each function, which specifies the elapsed time between the data -transmission- 
start and data -transmission-end events within the time profile. The effect of 

a difference between the input and output data transmission times is to change 
the time at which the data-transmission-end event occurs by the amount of 
that difference. The mechanics of these bookkeeping steps are straight- 
forward and identical for all processing functions and will not be treated 
further. 
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6.7 The input, output, and state-of-the-art limitations for each proces- 
sing function have been developed, assuming that the user will specify both the 
processing functions to be accomplished for each sensor and the order m which 
they are to be performed. As presently envisioned, the model will not consider 
the reasonableness of the processing requested or the order m which it is to be 
performed, with the exception that a check of the data type input to each pro- 
cessing function will be made to ensure that it is of the proper type (analog or 
digital) . 

6.8 The input and output for each processing function consist of three 
basic data items. These are' 

a. Information (analog data) or data (digital data) rate 

b. Ground information or data rate 

c. Time profile. 

In the ensuing development it is very important that the difference between the 
data rate (information rate) and the ground data rate (ground -information rate) be 
clearly understood. The data rate specifies the physical or actual input (or 
output) data rate to each processing function. The ground data rate specifies 
the minimum data rate required to transmit the information gathered' by the sensors 
in one complete sensor cycle time. In general the data rate can be related to the 
ground data rate by 

total ground datag 
data transmission time 

B 


data rate^ — 


total data^ 
total ground data. 


which becomes 

data rate 

B 


data ratej\ x data transmission timej\ 
ground data rate^ x cycle time 


-ground data rateg x cycle time 

data transmission time. 

L - c 


where 


data rate^ x data transmission iime^ 
ground data rate^ x cycle time 


accounts for the redundancy from resolution element to resolution element (or 
bit-to-bit), 

, cycle time 

and , ■ . ■ 

data transmission time D 

B 

L J 

accounts for the increase m data rate due to a reduction m the time alloted for 
transmission. 
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Removing the redundant terms this equation becomes 
data rate^ x data transmission timeA 


data rate^ = 
B 


ground data rate x data transmission time 
Pi B 


|ground data rateg 


6.9 The time profile specifies the time at which physical data transmission 
takes place and how that transmission is related to cycle and absolute time. 

This is accomplished by maintenance of a time profile containing an event list. 
T-he event list contains four possible types of events; these are* 

a. Cycle start 

' b. Data transmission start 

c. Data transmission end 

d. Cycle end. 

Signal Conditioning 

6.10 The signal conditioning function provides the capability to modify the 
analog information rate and the timing of the data transmission to the next pro- 
cessing function. Although it is not limited to such usage, the signal condition- 
ing function may be considered the link between the sensor and the output of the 
sensor system. 

6.11 Under the assumption that the analog information rate is related to the 
signal bandwidth by info rate = 2 x bandwidth, the increase or decrease m the 
output information rate of the signal conditioning function may be determined 
as a function of the input information rate and the ratio of output-to- input 
bandwidth. 

6.12 Since the signal conditioning function does not modify the information 
content of the input, the ground information rate is unchanged by the function, 
and the data transmission time (time required to transmit to the next processing 
function, that information collected by one cycle of the sensor) is increased or 
decreased by either the ratio of the input-to-output maximum bandwidth or by 
the ratio of actual mput-to-output bandwidth, for those cases m which the 
maximum capacity of the signal conditioning unit is not used. 

Analog/Digital Conversion 

6.13 The analog/digital conversion function provides the capability to 
convert the analog information rate to a digital data rate. Buffering is not 
considered to be inherently necessary within the A/D conversion function, since 
it may not be desirable in some applications, Therefore, it was not included. 

The output of the conversion may be transmitted to the next unit by either 
parallel or serial means, however, the number of bits transferred over a given 
period of time is, m general, unaffected by the means of transmission. Thus, 
the means of transmission has been ignored. 
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6.14 Under the assumptions stated previously, the output data rate of the 
A/D conversion function is given by 

data rate^ = (info rate^) log^ fo re y levels) . 

Although the A/D conversion does not modify the information content of the 
input, it is convenient for later comparison to convert the ground information 
rate to a ground daia rate, thus, 

grd data rate^ = (grd mfo rate^) log^ (grey levels) . 

Since buffering was not considered an integral part of the A/D conversion, 
no change m data transmission time occurs 

Multiplex or Commutate 

6.15 The multiplex or commutate function provides the capability of com- 
bining a number of separate inputs into a single output. The multiplex function 
operation is defined as follows If any input channel is active, the output is 
active at a constant data rate. Only m the case of no active input channel will 
the output be suppressed. Since each channel to the multiplex function has a 
separate and distinct time profile associated with it, the output of the multiplex 
function will be determined only after these time profiles have been merged and 
sufficient extra cycles are added to each to form a periodic function output. 

6.16 Based on the assumption of input- sampling rate synchronization, the 
output data rate is given by* 

data rate Q ~ (max data rate/channel) (number of channels). 

Since a number of periodic cycles for each input will be necessary to obtain a 
periodic output, the output ground data rate is 

n 

grd data rate^ = ^ (grd data rate^ ) (N^) 

i 

where grd data rate^. = ground data rate input for ith channel 

l 

N = number of cycles of the ith channel required 
1 for periodic output 

n = number of multiplexer channels . 

The data transmission times for the multiplex function will be determined 
logically from the periodic output time profile and will specify data transmission 
times for all times when the function output is active. 

Buffer 

6.17 The buffering function provides the capability to smooth the digital 
input by allowing modification to the data rate. Since, m general, it is 
desirable from an efficiency standpoint to match the output of the buffer with 
the maximum input data rate of the next function, the buffer function, within 
the specified state-of-the-art limitations, will match its output data rate to 
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that of the next function. The maximum output data rate for the buffer function 
is assumed to be equal to the maximum input data rate, since transfers into and 
out of storage may normally be expected to require equal time. 


6.18 Xhe output data rate of the buffer function is chosen as the maximum 
input data rate for the next function unless the state-of-the-art limitation for 
the buffer input data rate is exceeded, m which case the output is equal to the 
maximum buffer input data rate for the buffer function. Since the function does 
not modify the information content of the input, the ground data rate is unaffected 
by the function and the output data transmission time is increased or decreased 
in proportion to the ratio of input- to-output data rates. 

6.19 The storage required by the buffer function may be found by consider- 
ing the data received minus the data transmitted by the function during one cycle. 
Figure 8 depicts this situation for three cases: 

(1) The output data rate less than the input data rate 

(2) The output data rate equal to the input data rate 

(3J The output data rate greater than the input data rate. 

For case (1) the storage requirement is 

storage = (data rate^) (date transmission tirne^.) 

- (data rate^) (data transmission time^ - time delay). 

For cases (2) and (3) the storage requirement is 

storage = (data rate^) (time delay) . 

Case (2) is a special usage of a buffer function where the buffer is simply a 
time delay. Note that a trivial case exists (in addition to those listed above) , 
m which the time delay exceeds the input data transmission time. For this case 
the storage requirement is simply the total data received, i.e., the product of 
the input data rate and the input data transmission time. 



6.20 The overlap redundancy removal function provides the capability of 
removing duplicate data gathered by the sensor. The previous quarterly technical 
report discussed overlap redundancy by geometrical computation. Alternately, 
this data removal can be accomplished by comparison between data from two 
sensor cycles. For this case, a buffer will be included as an integral part of 
this function. Due to the required comparison, the buffer storage requirements 
can be expected to exceed the amount of storage required to store one cycle of 
data; therefore, the buffer is considered of sufficient size to allow matching of 
the output data rate with that of the next function within the stata-of-the-art 
limitations. 


6.21 The output data rate of the overlap redundancy removal function is, 
then, equal to the maximum input data rate for the next function, unless the 
state-of-the-art limitation for the overlap redundancy buffer is exceeded, m 
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which case the output is equal to the maximum input data rate for the overlap 
redundancy removal function. The amount of real information, the ground data 
rate over one cycle time, is reduced by the amount of overlap which has been 
removed. The ground overlap area in the following equation is the amount of 
overlap which is to be removed. This amount is not necessarily the full amount 
of image overlap. Recalling that the ground data rate has a fixed transmission 
time, the cycle time, the output ground data rate is given by 


grd data rate^ = grd data rate^ - 


grd ovrlp area 
grd res area 


log 2 (grey levels) 


x 


1 # 

cycle time 


The output data transmission time can be computed from the output data rate 
and the amount of data requiring transmission. The amount of data requiring 
transmission is given by the input data minus the removed redundant data; thus. 


data transmission time^ = 



(data rate T ) (data transmission time^) 

r\ 


- (eff 


data ratej x data transmission timej'. / grd ovrlp area 


grd data ratej x cycle time 


~i. 


grd res area 


x log 2 (grey levels) 


data rate 


where 


grd ovrlp area 
. grd res area 



accounts for the amount of overlap information removed, 


and 


'data ratei x data transmission timej 
i grd data rate^ x cycle time 


) 


accounts for the element-to-element redundancy. The redundancy removal 
efficiency is provided to compensate for both redundancy removal algorithm 
inefficiency and sensor-related problems such as sensor registration. 


6.22 Since the maximum input data rate for the overlap redundancy removal 
function is determined by the buffer, as is the' maximum output data rate (which 
is assumed to be equal to the maximum input data rate), it is necessary to re- 
quire a second state-of-the-art specification for the redundancy removal pro- 
cedures. Based on the assumption that the processing time required for removal 
of redundant data is a function of the amount of redundant data to be removed, a 
. maximum data removal rate provides this information. To ensure that the 
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redundant data can be removed without data buildup in the buffer and a resultant 
loss of data, the following inequality must hold: 



( data ratej x data xmit tirnej] (g rd ovr lP area'' 
grd data ratej x cycle time / y grd res area y 


(max data removal rate) 


l°g 2 


(grey levels) 


< cycle 
“time 


The overlap redundancy removal function time delay is also a function of the 
amount of data removed. Assuming that data removal always proceeds at the 
maximum rate, the additional time delay or processing time delay is> 


proc time delay = eff ! 


'data ratej x data xmit timej\ 7 grd ovrlp area 
, grd data ratej x cycle time./ \ grd res area y 


x 


] 


log2 (grey levels) 


max data removal rate 


It is conceivably possible, although unlikely, for the maximum data removal 
rate to exceed the maximum input data rate, m which case the maximum input 
data rate will be used m the processing time delay computation m place of the 
maximum data removal rate. 


6.23 The storage requirement for the overlap redundancy removal function 
can be computed on the same basis as the buffer function storage requirement, 
with the exception that a number of input data transmissions may occur prior 
to the occurrence of the case shown m Figure 8. Additional storage require- 
ments from these transmissions are computed by taking the product of the data 
transmission time and the data rate for each such occurrence and summing the 
products of all cases. 


Information Preserving Data Compression 

6.24 The information preserving data compression function provides the 
capability of removing the redundant data produced by processing of the sensor 
data collected. Since data compression does not, in general, produce a smooth 
output data rate, a buffer is included as an integral part of this processing 
function. To avoid placing an undue burden on the user, who must specify the 
state-of-the-art limitations for this processing function, the function of this 
buffer will be limited to smoothing of the output data rate, and the input and 
output data transmission times will be assumed to be equal. 


6.25 Recalling that the ground data rate is the minimum data rate achievable 
without loss of information and is smoothed over the cycle time, the output data 
rate is 


data rate^ - data rate^ - j^data rate^. - grd data rate^j |eff 

where grd data rate x = (grd data ratej) f-g— — ■ — irn ^ — 

\ I \data transmission timei 
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The compression efficiency is provided to compensate for the fact that the 
algorithm chosen may not completely remove the redundant information. Real 
information content is not changed; thus, the ground data rate is unchanged, and 
the data transmission time is unchanged (by assumption). The storage require- 
ment for the buffer is determined as it was for the buffer function. 

Entropy Reducing Data Compression 

6.26 The entropy reducing data compression function provides the capa- 
bility of selectively removing data to obtain significantly lower data rates 
while accepting the fact that some information may be lost as a consequence 

of the data reduction. As for the information preserving data compression func- 
tion, a buffer is included as an integral part of this function, and the output data 
rate is based on equal input and output data transmission times. Unfortunately, 
however, sufficient data concerning the algorithm to be used for data reduction 
and the amount of reduction possible due to the characteristics of the data col- 
lected are not available to the model as presently envisioned, for this reason, 
an anticipated reduction factor must be included with the state-of-the-art 
specifications for this function. 

6.27 The output data rate is 

, t data rater 

data rateQ = t . 

reduction factor 

The real information content is not changed, thus , the ground data rate is un- 
changed by this function, and the data transmission time is unchanged. 

6.28 The storage requirement for the integral buffer is determined as it 
was for the buffer function unless the anticipated buffer size as a percentage 
of the total input data is specified. If so, the buffer size is 

storage = (data rate^.) (data transmission tirne^.) (ant buf size %) . 

6.29 Table 1 presents, m tabular form, the input, output, state-of-the-art 
limitations data which must be supplied for each processing function. The 
methodology to be used m generation of the output data for each of the process- 
ing function has been discussed m the preceding paragraphs. 
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TABLE X 

PROCESSING FUNCTION SPECIFICATIONS 


Signal Conditioning 

Input: 

Information rate (analog) 

Ground information rate (analog) 

Time profile 

Output: 

Information rate (analog) 

Ground information rate (analog) 

Time profile 

State -of -the -art-data : 

Time delay (time from first input to first output) 
Maximum input bandwidth (Hz) 

Maximum output bandwidth ’■ (Hz) 

Ratio of output-to-input information rate (optional) 

Output computation 

Information rate: 

a . Info rate q 

b . Info rate^ 

Ground information rate 

Grd info rale Q = grd info rate-j- 


( max bandwidth q\ 
max bandwidth^. J in ^° rate j 

= (ratio of output-to -input info rate) info rate^ 

(optional— to be used if ratio is different from a. above) 


(resolution elements/sec) 
(resolution elements/sec) 
(sec) 

(resolution elements/sec) 
(resolution elements /sec) 
(sec) 


Time profile- 

a. Time profile = (time profiler) +■ time delay , -y 

O ^ i max bandwidthj 

b. Data transmission timeo — (data transmission time T ) — r — 3 — rrr — 

O' xy max bandwidthQ 

c. Data transmission timeQ = (data transmission timej) 

/ 1 ^ 


x 


— 1 

ratio of output-to -input info rat ej 


Analog/Dig ital Conversion 


Input 

Information rate (analog) 

Ground information rate (analog) 
Time profile 


(resolution elements/sec) 
(resolution elements/sec) 
(sec) 





TABLE 1 (Cont) 


Data rate (digital) (bits/sec) 

Ground data rate (digital) (bits/sec) 

Time profile (sec) 

Stare -of -the -art data : 

Time delay (time from first mpur to first output) 

Maximum input bandwidth (Hz) 

Maximum grey levels 
Desired grey levels (optional) 

Output computation. 

Data rale: 

a. Data rateQ = (info rate-jO log 2 (max grey levels) 

b. Data rateQ = (info rate^) log 2 (desired grey levels) 

(optional— to be used if desired grey 
levels are different from maximum) 

Ground data rate . 

a. Grd data rateQ = (grd info ratej) log2 (max grey levels) 

b. Grd data rateQ = (grd info ratej) log2 (desired grey levels) 

(optional) 

Time profile: 

a. Time profileQ= (time profilej) + time delay 

b. Data transmission timeQ = data transmission timej 


(Hz) 


Multiplex or Commutate 


Input (one to n inputs where n is the number of channels): 

Data rate (digital) (bits/sec) 

Ground data rate (digital) (bits/sec) 

Time profile (sec) 


Data rate (digital) 

Ground data rate (digital) 

Time profile 

State-of-the-art data : 

Time delay (time from first input to first output) 
Maximum input data rate/channel 
Output data rate (optional) 

Number of channels 


(bits/sec) 

(bits/sec) 

(sec) 
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TABLE 1 (Cont) 



Output computation. 

Data rate: 

a. Data rateQ = (max data rate/channel) (number of channels) 

b. Data rateQ = output data rate 

(optional — to be used if different from a. above) 

Ground data rate: 

n 

a . Grd data rate^ = £ (ground data ratej ) (N^ 

1 1 x 

where jN^ is the number of input data cycles 

Time profile: 

a. Additional cycles will be generated for each input until out- 
put is periodic . Output data rate is constant while any input 
is active . 

Buffer 

Input 

Data rate (digital) 

Ground data rate (digital) 

Time profile 

Output 

Data rate (digital) 

Ground data rate (digital) 

Time profile 

State-of-the-art data* 

Maximum time delay (time from first input to first output) 

Maximum input data rate 
Maximum storage available 
Time delay (optional) 

Output computation: 

Data rate 

a. Comparison with maximum input data rate for the next function 
will produce an output data rate of 

Data rateQ = max data ratej (next unit) 

if max data ratej (next unit) is less than the maximum input 
data rate of the buffer, or 

Data raieQ = max data ratej (buffer) 


(bits/sec) 

(bits/sec) 

(sec) 

(bits/sec) 

(bits/sec) 

(sec) 
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TABLE 1 (Cont) 


Ground data rate: 

a. Grd data rate^= grd data rate 
O I 

Time profile: 

a. Time profileQ = (time profile^.)-*- max time delay 

b. Time profile^ = (time profile^.) + time delay 
{optional — to be used if time delay other than 
maximum is desired) 

/ data rate 

c. Data transmission time,-, = (data transmission time ) 

° I ' data rate Q 


Overlap Redundancy Removal 

Input: 


Data rate (digital) 

(bits/sec) 

Ground data rate (digital) 

(bits/sec) 

Time profile 

(sec) 

Output* 


Data rate (digital) 

(bits/sec) 

Ground data rate (digital) 

(bits/sec) 

Time profile 

(sec) 


State-of-the-art data: 


Time delay (time from first input to first output with no removal of 
redundant data) 

Maximum data removal rate (function of the processing time required to 

locate and remove redundant data) 

Maximum input data rate 

Maximum storage available 

Redundancy removal efficiency 

Ground overlap area (supplied by overlap routine) 

Ground resolution area (from sensor routine) 

Grey levels (from sensor or A/D conversion) 

Output computation: 

Data rate: 

Since a buffer is an integral part of this processing function, the output 
data rate will be matched to the input data rate of the next unit. 
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r- 


a. Data rate^= data rate^ (next unit) 

if data rate next unit) is less than maximum output data rate, or 

b. Data rateQ- maximum data ratej 
Ground data rate 

~ , . , , / ard ovrlp area 

a. Grd data rate,= grd data rate., - — ; 

O I V 9 r d area 


log^ (grey levels) 


x 


cycle time 

Time profile 

a. Time profile^ = (time profile^.) + time delay + proc time delay 


(Eff) y 


where proc time delays 

/ data rate^. x data transmission time^ / grd ovrlp area \ 


grd data rate^ x cycle time 


( grd ovrlp area \ 

— C /log. (grey levels) 

grd res area J 2 


(max data rem rate) 


b. Data transmission time^ = 


(data rate^} (data transmission time^.) - (Eff) 


/ qrd ovrlp area) , . , . . ( 1 \ 

x log^ (grey levels) i" / 

\ grd res area J 2 J\data rat e^y 


Information Preserving Data Compression 

Input* 


Data rate (digital) 

(bits/sec) 

Ground data rate (digital) 

(bits/s ec) 

Time profile 

(sec) 

Output: 


Data rate (digital) 

(bits/s ec) 

Ground data rate (digital) 

(bits/sec) 

Time profile 

(sec) 

State-of-the-art data: 


Time delay (time from first input to first output) 

Maximum input data rate 
Maximum storage available 
Compression efficiency 



data rate^ x data transmission tirne^ 
v grd data rate^. x cycle time 
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Output Computation: 

Data rate: 

a. data rate Q - data rate^. - (data rate^. - grd data rate^) (eff) 


where 

grd data rate = fgrd data rate 


' / cycle time 

I/ldata transmission time. 


Ground data rate: 

a. Grd data rateQ= grd data rate^ 

Time profile 

a . Time profile Q = (time profile^.) + time delay 

b. Data transmission timeQ= data transmission time^ 


Entropy Reducing Data Compression 


Input 


Data rate (digital) 

(bits /sec) 

Ground data rate (digital) 

(bits /sec) 

Time profile 

(sec) 

tput 

Data rate (digital) 

(bits /sec) 

Ground data rate (digital) 

(bits /sec) 

Time profile 

(sec) 


State-of-the-art data 

Time delay (time from first input to first output) 

Maximum input data rate 

Anticipated buffer size {% of input data) (optional) 

Maximum storage available 

Reduction factor (input to output reduction) 

Output computation: 

Data rate 

a. Data rate^^ data rate^/reduction factor 
Ground data rate 

a . Grd data rateQ = grd data rate^ 

Time profile 

a . Time profile^ = (time profile^) + time delay 

b. Data transmission time = data transmission time T 

O I 
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APPENDIX A 

TIME-IN-VISUAL- VIEW MODEL 
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FIGURE A . 1 . (Coni) 
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FIGURE A. 1. (Cont) 
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B = BETA (LONGTB, 

LATTB,LONGES(lES), 

LATES(IES)) 



FIGURE A. 1. (Cont) 
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FIGURE A. 1. (Cont) 
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FIGURE A. 1. (Cont) 
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FIGURE A . 2 . LONGITUDE FUNCTION ROUTINE 

A-12 
















FIGURE A. 2. (Cont) 
A-13 










FIGURE A. 2. (Cont) 
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FIGURE A .3 . BETA FUNCTION ROUTINE 
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TABLE A. 1 

INPUT AND OUTPUT FOR TIME -IN -VISUAL -VIEW ROUTINES 


Titne-In-Visual View Routine 

Input 



1 Orbital Period (minutes) 

— 

PO 

2 Earth Station Latitude (radians) [Array] 

- 

LATES 

3 Earth Station Longitude (radians) [Array] 

- 

LONGES 

4 Longitude of First Ascending Node (radians) 

- 

LONGO 

5 Orbit Inclination Angle (radians) 

- 

AI 

6 Earth Station Minimum Look Angle (radians) [Array] 

- 

ALOOKM 

7 Satellite Altitude (nmi) 

- 

H 

8 Number of Orbits in Schedule Time Period 

— 

N 

9 Number of Earth Stations 

- 

M 

Output 



1 Satellite Tune of Arrival at Earth Station (minutes) 

— 

TARVL 

2 Satellite Time m Visual View (minutes) 

- 

TIV 

3 Earth Station Identification 


IES 

Longitude Function Routine 

Input 



1 Satellite Rotation Angle (radians) 

— 

X 

2 Orbit Inclination Angle (radians) 

- 

AI 

3 Orbital Period (minutes) 

- 

PO 

4 Longitude of First Ascending Node (radians) 

— 

LONGO 

Output 



1 Satellite Longitude (radians) 

- 

LONGF 

Beta Function Routine 

Input 



1 Satellite Longitude (radians) 

— 

SLN 

2 Satellite Latitude (radians) 

- 

SLT 

3 Earth Station Longitude (radians) 

- 

ESLN 

4 Earth Station Latitude (radians) 


ESLT 

Output 



1 Distance Between SuborbitaLPomt and Earth Station (radians) 

- 

BETA 
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APPENDIX B 
EARTH MODEL 


B-l 



B-2 


SUB POINT T RACK 



FIGURE B . 1 . EARTH MODEL SLIDE RULE 


APPENDIX C 

SPACECRAFT DATA PROCESSING MODEL— PROCESSING 
FUNCTION FLOW CHARTS 
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FIGURE C . 1 . SIGNAL CONDITIONING ROUTINE 
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FIGURE C . 1 (Cont) 



C-3 





















FIGURE C. 2 (Cont) 
C-6 














FIGURE C .3 . MULTIPLEX OR COMMUTATE ROUTINE 


0-7 











N/ 


S 

Cycle = Cyc Tj{1) 
Njd) = 1 
K= 1 


Note Cyc T j is an array containing the 

input cycle times of each multiplier 
mput 

Nj is an array containing the number of 
cycles of each input for a periodic 
output 



ISW= 1 
A = Cyc Tj(K) 
B = Cycle 


Yes 


Is Cyc Tj (K) > 
Cycle 7 







Multiplexer 
Input Cannot be 
Made Periodic 




FIGURE C. 3 (Cont) 
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FIGURE C. 3 (Cont) 
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FIGURE C. 3 (Cont) 
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BUFFER 


Availablej = May Data Ratej 
— Data Ratej 


Data RateQ = 
Max Data Ratej 
(Next Unit) 







T 


Grd Data Rate Q = 
Grd Data Ratej 


1 

Data Xmit Time 0 = 
(Data Xmit Timej) 

1 J Data Ratej \ 

\ Data Rate 0 j 


I 



FIGURE C . 4 . BUFFER ROUTINE 
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FIGURE C. 4 (Cont) 
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FIGURE C. 4 (Cont) 
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FIGURE C .5 . OVERLAP REDUNDANCY REMOVAL ROUTINE 
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X 


Residual = Time Delay + Proc 


Residual = Tune Delay 

Time Delay — (N)(Cyele Time) 


t Proe Time Delay 



FIGURE C . 5 (Cont) 
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FIGURE C. 5 (Cont) 
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FIGURE C. 5 (Cont) 
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FIGURE C 6 . INFORMATION PRESERVING DATA COMPRESSION ROUTINE 
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FIGURE C . 6 (Cont) 
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J = 0 
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J = J + 1 


Time Profile 0 (J) = 



FIGURE C. 6 (Cont) 
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FIGURE C. 7 ENTROPY REDUCING DATA COMPRESSION ROUTINE 
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FIGURE C. 7 (Cont) 
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FIGURE C. 7 (Cont) 
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